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2^ ■ Abstract 

Qj ', The Relativistic Hartree Bogoliubov (RHB) theory is apphed in the descrip- 

Q I tion of ground-state properties of Ni and Sn isotopes. The NL3 parameter set 

Q>^ ' is used for the effective mean-field Lagrangian, and pairing correlations are 

described by the pairing part of the finite range Gogny interaction DIS. Fully 

►-j^ . self-consistent RHB solutions are calculated for the Ni (28 < N < 50) and Sn 

Tlj" i (50 < N < 82) isotopes. Binding energies, neutron separation energies, and 



proton and neutron rms radii are compared with experimental data. The 

^SJ ' model predicts a reduction of the spin-orbit potential with the increase of 

^_^ the number of neutrons. The resulting energy splittings between spin-orbit 

Q^ . partners are discussed, as well as pairing properties calculated with the finite 

range effective interaction in the pp channel. 
^^ ■ 

g ■ PACS numbers: 21.10.D, 21.10.F, 21.60.J, 21.30, 27.50, 27.60 

> 

X 

H ■ I. INTRODUCTION AND OUTLINE OF THE MODEL 

Relativistic models of nuclear structure provide a consistent framework in which the 
nuclear many-body system is described in terms of interacting baryons and mesons. Detailed 
properties of finite nuclei along the /3-stability line have been very successfully described 
in the framework of relativistic mean field models (for a recent review see Ref. []I|]). In 
addition to the single-nucleon mean-field, pairing correlations have to be taken into account 
for a quantitative description of ground-state properties of open-shell nuclei. For nuclei 
close to the /3-stability line, pairing has been included in the relativistic mean-field model 
in the form of a simple BCS approximation p|, with monopole pairing force adjusted to 
the experimental odd-even mass differences. Pairing correlations are also crucial for the 
description of deformations in heavy nuclei. However, as we move away from the valley 
of /9-stable nuclei, the ground-state properties calculated within the BCS approximation 
become unreliable. In particular properties that crucially depend on the spatial extensions 
of nucleon densities, as for example nuclear radii. The reason is that the BCS scheme does 



not provide a correct description of the scattering of nucleonic pairs from bound states to 
the positive energy particle continuum [^^ , and it leads to unbound systems. The solution, 
of course, is a unified description of mean-field and pairing correlations, as for example 
in the framework of the Hartree-Fock-Bogoliubov (HFB) theory. Within a non-relativistic 
approach to the nuclear many-body problem, HFB theory has a long and successful history 
of applications. In particular, HFB theory in coordinate space [H has been used to describe 
properties of exotic nuclei with extreme isospin values, both on the neutron rich side [Q, 
and proton drip-line nuclei 0. 

The relativistic extension of the HFB theory was originally derived by Kucharek and 
Ring 1^. Starting from the Lagrangian of quantum hadrodynamics, they have been able 
to show that the pairing correlations result from the one-meson exchange {a-, uo- and p- 
mesons). The Relativistic Hartree-Bogoliubov (RHB) model developed in Ref. is based 
on the Hartree approximation for the self-consistent mean field. However, if for the one- 
meson exchange pairing interaction in RHB one uses the coupling constants from standard 
parameter sets of the RMF model, the resulting pairing correlations are much too strong. 
The repulsion produced by the exchange of vector mesons at short distances results in a 
pairing gap at the Fermi surface that is by a factor three too large. On the other hand, it 
has been argued in many applications of the Hartree-Fock-Bogoliubov theory that there is no 
real reason for using the same effective forces in both the particle-hole and particle-particle 
channel. In a first-order approximation, the effective interaction contained in the mean-field 
r is a G-matrix, the sum over all ladder diagrams. The effective force in the pp channel, i.e. 
in the pairing potential A, should be the K matrix, the sum of all diagrams irreducible in 
j9j9-direction. However, very little is known about this matrix in the relativistic framework 
and only phenomenological effective forces have been used in the pp channel of RHB. In the 
RHB calculations of Ref. pairing correlations have been described by a two-body force 
of finite range of Gogny type |§], 

yPP(l,2) = ^ e-«"^-"^)/^')'(iyi + £iP"-ifiP"-MiP"P"), (1) 

i=l,2 

with the parameters /Xj, VFj, Pj, Hi and Mj {i = 1, 2). For the DIS ||^ parameter set of the 
interaction the model was applied in the study of several isotope chains of spherical Pb, Sn 
and Zr nuclei. The pairing interaction is a sum of two Gaussians with finite range and prop- 
erly chosen spin and isospin dependence. The Gogny force has been very carefully adjusted 
to the pairing properties of finite nuclei all over the periodic table. Its basic advantage is 
the finite range, which automatically guarantees a proper cut-off in momentum space. In 
Refs. P,pi]| we have used RHB in coordinate space with the DIS Gogny interaction to study 



properties of light nuclear systems (C, N, O, F, Ne, Na, Mg) with large neutron excess. 
Self-consistent solutions were calculated for the ground-states of a number of neutron-rich 
nuclei. Predictions were obtained for the location of the neutron drip-line, reduction of the 
spin-orbit interaction, rms radii, changes of surface properties, formation of neutron skin 
and neutron halo. The effective interactions in the pp channel of the HFB theory were 
recently discussed in Ref. 0]. In particular, the role of finite range and the importance of 
density dependence were analyzed. Finite range forces have the obvious advantage of the 
automatic cut-off of high momentum components. On the other hand computations with 
zero-range forces are much simpler, although an artificial energy cut-off has to be included. 



The density dependence of interactions in the pairing channel influences the spatial distri- 
butions of pairing densities and flelds. By using interactions without density dependence, 
strong pairing flelds are produced in the volume of the nucleus. There is experimental evi- 
dence that pairing is a surface effect, and by including a density dependent component the 
pairing fleld can be made surface peaked. A fully self-consistent RHB model in coordinate 
space, with a density dependent interaction of zero-range (delta force), has been used to 
describe the two-neutron halo in ^^Li [|TT|. 
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In going away from the valley of /3-stable nuclei, the main problem of nuclear struc- 
ture models becomes the extrapolation of effective forces to nuclei with extreme isospin 
values. We have already applied RHB with the Gogny force in the pairing channel to light 
neutron-rich nuclei. However, many of these nuclei are still not accessible in experiments 
and therefore many of our results could not be compared with empirical data. In order to 
make predictions for medium-heavy nuclei at the neutron drip-line, we have to test available 
effective interactions in detailed calculations of properties of neutron rich nuclei for which 
a comparison with experimental data is possible. In the present article we consider two 
sets of isotopes: Ni (28 < N < 50) and Sn (50 < A^ < 82). Such an analysis will test 
the predictions of effective forces in both the ph and pp channel over two major neutron 
shells. In addition to the effects of the pairing interaction, we are particularly interested 
in the behavior of the spin-orbit term of the effective potential as function of the neutron 



number. For light neutron rich nuclei, in Ref. |T2| we have shown that the magnitude of 
the spin-orbit potential is considerably reduced at the drip-line, resulting in much smaller 
energy splittings between spin-orbit partners. For the Ne isotopes this reduction was found 
to be around forty percent. Since it seems that at present only relativistic models include 
the correct isospin dependence of the spin-orbit term in the mean-fleld potential, it would be 
important to study in medium-heavy nuclei the predicted spacings of single-neutron levels 
close to the Fermi surface. 



The details of the RHB model are given in Refs. 110,13|. A short outline of essential 
features follows. The ground state of a nucleus |$ > is described as vacuum with re- 
spect to independent quasi-particle operators, which are deflned by a unitary Bogoliubov 
transformation of the single-nucleon creation and annihilation operators. The generalized 
single-nucleon Hamiltonian contains two average potentials: the self-consistent mean-fleld 
r which encloses all the long range ph correlations, and a pairing fleld A which sums up 
the pp correlations. The expectation value of the nuclear Hamiltonian < $|i5f|$ > is a 
function of the hermitian density matrix p, and the antisymmetric pairing tensor k. The 
variation of the energy functional with respect to p and k, produces the single quasi-particle 



Hartree-Fock-Bogoliubov equations ||I^. The relativistic extension of the HFB theory was 
introduced in Ref. [Q. In the Hartree approximation for the self-consistent mean fleld, the 
Relativistic Hartree-Bogoliubov (RHB) equations read 
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where ha is the single-nucleon Dirac Hamihonian @) and m is the nucleon mass. The 
chemical potential A has to be determined by the particle number subsidiary condition, in 
order that the expectation value of the particle number operator in the ground state equals 
the number of nucleons. The column vectors denote the quasi-particle wave functions, and 
Ek are the quasi-particle energies. The Dirac Hamiltonian 

{\ — T-) 
ho = -ia-V + f3{m + g,(r)+g^ij'> + gpT3pl + e ^ ^ ^M ° (4) 

contains the mean-field potentials of the isoscalar scalar cr-meson, the isoscalar vector uu- 
meson and the isovector vector p-meson. A'^ is the electrostatic potential. The RHB equa- 
tions have to be solved self-consistently, with potentials determined in the mean-field ap- 
proximation from solutions of Klein-Gordon equations 

hA + ml] a(r) = -g^ p,(r) - g2 a\r) - g-, a^r) (5) 

[-A + ml]u\r) = -g^p,ir) (6) 

[-A + ml]p\r) = -g,p,{r) (7) 

-AA\r) = ep,ir). (8) 

for the sigma meson, omega meson, rho meson and photon field, respectively. The spatial 
components uj, p^ and A vanish due to time reversal symmetry. The equation for the sigma 



meson contains the non-linear a self-interaction terms |15[. Because of charge conservation, 
only the 3-component of the isovector rho meson contributes. The source terms in equations 
(I) to (H) are sums of bilinear products of baryon amplitudes 

Ps(r) = E ^iW7V,(r), (9) 

Ek>Q 

P.(r) = E ^iW^.W, (10) 

P3(r) = E V^{v)nV,{Tl (11) 

Ek>0 

pUr) = E yHr)^V,{r). (12) 

where the sums run over all positive energy states. The pairing field A in (^ is defined 

A,,(r, r') = ^ E ^afecd(r, r')«:,,(r, r'). (13) 

where a, 6, c, d denote quantum numbers that specify the single-nucleon states. Vahcdiji r') 
are matrix elements of a general two-body pairing interaction, and the pairing tensor is 
defined 

«:,,(r,r')= E f/4(r)l^*(r'). (14) 

Sfe>0 



The eigensolutions of Eq. (^ form a set of orthogonal and normahzed single quasi-particle 
states. The corresponding eigenvalues are the single quasi-particle energies. The self- 
consistent iteration procedure is performed in the basis of quasi-particle states. The resulting 
quasi-particle eigenspectrum is then transformed into the canonical basis of single-particle 
states, in which the RHB ground-state takes the BCS form. The transformation determines 
the energies and occupation probabilities of the canonical states. 

For nuclear systems with spherical symmetry the fields cr{r), uj^{r), p^{r), and A^[r) 
depend only on the radial coordinate r. The nucleon spinors Uk (V^) in (|^) are characterized 
by the angular momentum j, its ^-projection m, parity vr and the isospin ^3 = ±| for neutron 
and proton. The two Dirac spinors Uk{r) and Vfc(r) are defined 

U.iVMr.s.t,) = Hr'^dr^^^^'^^XrM- (15) 

\^JUiV){r)^ij^l^rn{^,ip,S) J 

g{r) and /(r) are radial amplitudes, Xr is the isospin function, and VLji^ is the tensor product 
of the orbital and spin functions 

ms,mi ^ 

The two-component functions 

*„W:^(-«) and *.«:^(-«). (17) 

are solutions of the Dirac-Hartree-Bogoliubov equations 
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The self-consistent solution of the Dirac-Hartree-Bogoliubov integro-differential eigenvalue 
equations and Klein-Gordon equations for the meson fields determines the nuclear ground 
state. In Refs. [y,|l^,|l^,|l^] we have used Finite Element Methods in the coordinate space 
discretization of the coupled system of equations. In order to correctly describe structure 
phenomena in exotic nuclei with extreme isospin values, as for example regions of neutron 
halos with very diffuse neutron densities, the RHB equations have to be solved in coordinate 
space. However, for the spherical nuclear systems that we consider in the present article this 
is not necessary |0 . As it was done in Ref . J^ , the Dirac-Hartree-Bogoliubov equations and 
the equations for the meson fields are solved by expanding the nucleon spinors Uk{r) and 
Vfc(r), and the meson fields in a basis of spherical harmonic oscillators for A^ = 20 oscillator 



shells 0. For odd-A isotopes we also include the blocking procedure ||1J] for neutron levels. 
Strictly speaking, for these nuclei also the assumption about currents is no longer valid, and 
in general one should include the fields a;, p^ and A. However, the spatial components of 
the vector fields do not significantly contribute to the ground state properties investigated 
in the present work , and therefore they have not been included in our calculations. 



II. GROUND-STATES OF NI AND SN ISOTOPES 

The details of ground-state properties of Ni (28 < A^ < 50) and Sn (50 < iV < 82) 
isotopes will depend on the choice of the effective mean-field Lagrangian in the ph channel, 
as well as on the effective pairing interaction. Several parameter sets of the mean-field 
Lagrangian have been derived that provide a satisfactory description of nuclear properties 
along the /3-stability line. In particular, NLl |T^, NL3 |17|, and NL-SH [18|. The effective 



forces NLl and NL-SH have been used in many analyses to calculate properties of nuclear 
matter and of finite nuclei, and have become standard parameterizations for relativistic 
mean-field calculations. The effective interaction NLl was also used in the RHB+Gogny 
calculations of Ref. [0. More recently the parameter set NL3 has been derived [1^ by 



fitting ground state properties of a large number of spherical nuclei. Properties calculated 
with the NL3 effective interaction are found to be in very good agreement with experimental 
data for nuclei at and away from the line of /3-stability. In Ref. |jT9| it has been shown that 
constrained RMF calculations with the NL3 effective force reproduce the excitation energies 
of superdeformed minima relative to the ground-state in ^^^Hg and ^^^Pb. In the same 
work the NL3 interaction was also used for calculations of binding energies and deformation 
parameters of rare-earth nuclei. We have used the NL3 parameter set in our study of 
light neutron- rich nuclei in Refs. P JTIl| , p!^ . For Ni and Sn the objective is to study how 
well does the NL3 effective force describe systems with very different number of neutrons, 
without going to nuclei with extreme isospin values. Only if it were shown that the isospin 
dependence of NL3 is correct, this interaction could be used to make reliable predictions 
about medium-heavy drip-line nuclei. This could be especially important for Ni, since 
there is hope that drip-line isotopes might become accessible in experiments. Pairing being 
essentially a non-relativistic effect, we use the Gogny interaction in the pp channel with the 
parameter set DIS [Q. Results obtained with this effective force might also indicate the 
path one should take in deriving a fully relativistic theory of pairing, consistent with the 
modern mean-field Lagrangians. 

In Figs. |l| and ^ we display the one- and two-neutron separation energies 

Sn{Z, N) = B^{Z, N) - B^{Z, N-1) (19) 

S^niZ, N) = B^{Z, N) - B^{Z, N-2) (20) 

for Ni (24 < N < 50) and Sn (50 < N < 88) isotopes, respectively. The values that corre- 
spond to the self-consistent RHB ground-states are compared with experimental data and 



extrapolated values from Ref. |^ . The theoretical values reproduce in detail the experimen- 
tal separation energies. The model describes not only the empirical values within one major 
neutron shell, but it also reproduces the transitions between major shells. The results are 
excellent for the region beyond the shell closure at A^ = 82 in Sn. The agreement with ex- 
perimental data is somewhat worse for neutrons in the lf7/2 orbital in Ni isotopes, although 
the general trend is reproduced. However, these nuclei (24 < N < 28) really belong to the 
proton-rich side and are not a primary objective of the present study. The total binding 
energies for Ni and Sn isotopes are compared with experimental values in Fig. 0. Except 
for the region around ^°Ni and for ^°''~^°^Sn, the absolute differences between the calculated 
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and experimental masses are less than 2 MeV. For Ni the model predicts weaker binding for 
A^ < 40. Compared to experimental values, the theoretical binding energies are ~ 1 MeV 
larger for neutrons in the lg9/2 orbital (40 < N < 50). For Sn isotopes the results of model 
calculation in general display an overbinding. The differences are larger for ^''°~^''^Sn, but 
we expect that for these nuclei additional correlations should be taken into account in order 
to get a better agreement with experimental data. In particular, proton-neutron pairing 
could influence the masses in this region. Proton-neutron short-range correlations are not 
included in our model. 

In Fig. ^ we show the self-consistent ground-state neutron densities for the Sn and Ni 
nuclei. The density profiles display shell effects in the bulk and a gradual increase of neutron 
radii. In the insert of Fig. ^ we include the corresponding values for the surface thickness 
and diffuseness parameter. The surface thickness t is defined to be the change in radius 
required to reduce p{r)/ p^ from 0.9 to 0.1 (po is the maximal value of the neutron density; 
because of shell effects we could not use for po the density in the center of the nucleus). The 
diffuseness parameter a is determined by fitting the neutron density profiles to the Fermi 
distribution 

p(r) = Po (l + ea;p(^^)) \ (21) 

where Rq is the half-density radius. By adding more units of isospin the value of the neutron 
surface thickness increases and the surface becomes more diffuse. The increase in t and a is 
more uniform in Sn, and both parameters increase approximately forty percent from ^°°Sn 
to ^^^Sn. A somewhat smaller increase in the surface thickness is observed for Ni isotopes. 
The diffuseness parameter for Ni is essentially a step function: a ~ 0.45 fm for A^ < 40 
and a ~ 0.50 fm for neutrons in the lg9/2 orbital. We will show that the observed changes 
in surface properties result from the reduction of the spin-orbit term in the effective single- 
nucleon potential. 

In Figs. ^ and ^ we display the neutron and proton rms radii for Ni and Sn isotopes, 
respectively. The calculated values are compared with experimental neutron radii from Ref. 
pl| , and with data for proton radii from Ref. [^]. In the lowest panels we also compare 



the differences r„ — Vp. We find an excellent agreement between experimental data and 
values calculated with the NL3 effective force with the DIS Gogny interaction in the pairing 
channel. The model predicts a uniform increase of rms radii with the number of neutrons. 
The neutron skin r„ — r^ increases to approximately 0.4 fm at the closed shells A^ = 50 for 

Ni, and A^ = 82 for Sn. 



In Ref. [|T2| we have shown that in the framework of the relativistic mean-field model the 
magnitude of the spin-orbit term in the effective single nucleon potential is greatly reduced 
for light neutron rich nuclei. With the increase of the number of neutrons the effective spin- 
orbit interaction becomes weaker and this results in a reduction of the energy splittings for 
spin-orbit partners. The reduction in the surface region was found to be as large as ~ 40% 
for Ne isotopes at the drip-line. The spin-orbit potential originates from the addition of two 
large fields: the field of the vector mesons (short range repulsion), and the scalar field of 
the sigma meson (intermediate attraction). In the first order approximation, and assuming 
spherical symmetry, the spin orbit term can be written as 



Vs.o. = -|-V5s(r), (22) 

r or 



where Vis is the spin-orbit potential [p^ 
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Vis = {V~S). (23) 

V and S denote the repulsive vector and the attractive scalar potentials, respectively, me// 
is the effective mass 

m,ff = m-\{y-S). (24) 

Using the vector and scalar potentials from the NL3 self-consistent ground-state solutions, 
we have computed from (^2]) - (|2^) the spin-orbit terms for the Ni and Sn isotopes. They 
are displayed in Figs. |^ and |^ as function of the radial distance from the center of the 
nucleus. The magnitude of the spin-orbit term Vs.o. decreases as we add more neutrons, i.e. 
more units of isospin. If we compare ^^Ni with ^^Ni, in Fig. |^, the reduction is ~ 35% in 
the surface region. This implies a significant weakening of the spin-orbit interaction. The 
minimum of Vs.o. is also shifted outwards, and this reflects the larger spatial extension of 
the scalar and vector densities, which become very diffuse on the surface. The reduction 
of the spin-orbit term seems to be less pronounced in the Sn isotopes (Fig. ^, and this 
indicates that the weakening of the spin-orbit interaction might be not that important in 
heavy nuclei. The effect is reflected in the calculated spin-orbit splittings of the neutron 
levels in the canonical basis 

^Eis = En,lJ=l-l/2 — En,lJ=l+l/2, (25) 

In Fig. ^ we display the energy splittings of spin-orbit neutron partners for Ni and Sn, 
respectively. The calculated spacings are shown as function of the neutron number. We 
only include the spin-orbit doublets for which one of the partners is an intruder orbital in 
a major shell. These doublets display the largest energy splittings. We notice in Fig. |^ 
that the spacing between spin-orbit partners decreases with neutron number. The effect is 
stronger in Ni than in Sn, and quantitatively the observed decrease is consistent with the 
gradual weakening of the spin-orbit term shown in Figs. |^ and ^. 

In order to illustrate the properties of the interaction in the pp channel, in Figs. |TD| and ^ 
we plot the average values of the neutron canonical pairing gaps Anij as functions of canonical 
single-particle energies. The gaps are displayed for canonical states that correspond to the 
self-consistent ground states of ™Ni and ^^°Sn, respectively. Anij are the diagonal matrix 
elements of the pairing part of the RHB single-nucleon Hamiltonian in the canonical basis. 
Although not completely equivalent, Anij corresponds to the pairing gap in BCS theory. 
A very detailed discussion of HFB equations in the canonical basis can be found in Ref. 
Q. For ^^°Sn (Fig. 0) we essentially reproduce the results calculated in Ref. [^] with the 
Gogny HFB+DIS model. The calculated pairing gaps are almost identical, although there 
is a difference in the canonical single-neutron energies due to different effective forces used 
in the ph channel. While we employ the NL3 effective meson-exchange interaction in the ph 
channel, in the calculation of Ref. [^ the Gogny effective force has been used both in the ph 
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and pp channel. The pairing gaps display a uniform decrease with single-neutron energies. 
This is related to the volume character of the Gogny interaction in the pairing channel. The 
values of the Anij for the Si/2 orbitals are slightly larger, and the average value at the Fermi 
surface is around 1.5 MeV. It should be noted that for the HFB+SkP calculations of Ref. 
the values of the pairing gaps are strongly peaked at the Fermi surface. In addition, we 
have calculated the average canonical pairing gaps for ™Ni (Fig. |10|). As compared with 
^^°Sn, the Anij display a very similar behavior for deep-hole states (although for lsi/2 the 
pairing gap is considerably smaller), but they decrease more slowly for canonical states in 
the single-neutron continuum. At the Fermi surface the average values are between 1 and 
1.5 MeV, somewhat smaller than in ^^''Sn. 



Finally in Fig. |T2] we compare the averages of the neutron pairing gaps for occupied 
canonical states 

< A^ >= ^"'^-^"f-^ (26) 

2-inlj ^nlj 

where v'^ij are the occupation probabilities. The quantities < A^ > are plotted as functions 
of the neutron number for Ni and Sn isotopes. Solid lines correspond to even, and dashed 
lines to odd-neutron isotopes. < A^v > provides an excellent quantitative measure of pairing 
correlations. The quasi-parabolic functional dependence on the number of neutrons reflects 
the increase of pairing correlations toward the middle of a major shell. Both for Ni and 
Sn the values of < A^- > are above 2 MeV in the middle of the corresponding shells. As 
one would expect, the values for odd- neutron isotopes are somewhat lower. While in Sn the 
average pairing gaps do not provide any indication of additional shell effects, in Ni they very 
clearly indicate the shell subclosure at A^ = 40. It is also interesting to note that pairing 
correlations are stronger in the subshell 28 < A^ < 40, than for neutrons in the lg9/2 orbital. 
In conclusion, in the present work we have performed a detailed analysis of ground-state 
properties of Ni (28 < N < 50) and Sn (50 < N < 82) nuclei in the framework of the 
Relativistic Hartree Bogoliubov theory. The NL3 parameter set has been used for the effec- 
tive mean-field Lagrangian in the ph channel, and pairing correlations have been described 
by the finite range Gogny interaction DIS. In a comparison with available experimental 
data, we have shown that the NL3 -|- Gogny DIS effective interaction provides an excellent 
description of binding energies, neutron separation energies, and proton and neutron rms 
radii. The results indicate that this choice of model parameters might also be valid for nuclei 
with more extreme isospin values, i.e. medium-heavy nuclei at the drip-lines. We have also 
discussed the predicted reduction of the effective spin-orbit potential with the increase of the 
number of neutrons, as well as the resulting energy splittings between spin-orbit partners 
and modifications of surface properties. Pairing properties calculated with the finite range 
effective interaction in the pp channel have been carefully analyzed. These results are partic- 
ularly important since one of the main objectives of RHB theory should be a fully relativistic 
description of pairing correlations, consistent with modern mean-field Lagrangians. 
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FIGURES 

FIG. 1. One and two-neutron separation energies for Ni isotopes calculated in the RHB model 
and compared with experimental data. 

FIG. 2. Comparison between RHB model and experimental one and two-neutron separation 
energies for Sn isotopes. 

FIG. 3. Differences between RHB model and experimental binding energies for Ni and Sn 
isotopes. 

FIG. 4. Self-consistent RHB single-neutron density distributions for Sn (50 < N < 82) and 
Ni (28 < N < 50) nuclei, calculated with the NL3 effective interaction. 

FIG. 5. Calculated neutron and proton rms radii for Ni isotopes compared with experimental 
data. 



FIG. 6. Calculated and experimental neutron and proton rms radii for Sn isotopes. 

FIG. 7. Radial dependence of the spin-orbit term of the potential in self-consistent solutions 
for the ground-states of Ni (28 < N < 50) nuclei. 

FIG. 8. Radial dependence of the spin-orbit term of the potential in self-consistent solutions 
for the ground-states of Sn (50 < N < 82) nuclei. 

FIG. 9. Energy splittings between spin-orbit partners for neutron levels in Ni and Sn isotopes, 
as functions of neutron number. 

FIG. 10. Average values of the neutron canonical pairing gaps as functions of canonical 
single-particle energies for states that correspond to the self-consistent ground state of ^''Ni. The 
NL3 parametrization has been used for the mean-field Lagrangian, and the parameter set DIS for 
the pairing interaction. 

FIG. 11. Same as in the previous figure, but for the canonical states that correspond to the 
ground state of ^^'^Sn. 

FIG. 12. Average neutron pairing gaps < An > for the Ni and Sn isotopes, as functions of 
neutron number. 
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